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Abstract 

Non-perturbative effect of tlie formation of a chiral symmetry breaking condensate 
{ipip) and of a dynamically generated fermion mass in QED in the presence of an exter- 
I nal magnetic field is considered. The dynamical mass of a fermion (energy gap in the 

fermion spectrum) is shown to depend essentially nonanalytically on the renormalized 



o 
o 



ps) ■ coupling constant a in a strong magnetic field. Possible applications of this effect are 



o 
o 



discussed. 



The dynamics of fermions in a strong external magnetic field has been attracting much 
attention during last years. Perhaps, the brightest example has been the discovery and 
theoretical explanation of the fractional Hall effect leading to the 1998 Nobel Prize award 
j> ', (see Nobel lectures by Laughlin, Stormer and Tsui in [|l|) for "discovery that electrons acting 
O ■ together in strong magnetic fields can form new types of "particles", with charges that are 
^ fractions of electron charges ", as is said in press release of the Royal Swedish Academy of 
Sciencies. Thus, strong magnetic fields can drastically affect the ground state of a system 
leading to new types of excitations. 
O ■ In this talk, I will describe one more phenomenon in an external magnetic field: dynamical 
breaking of chiral symmetry induced by such a field, hence the name magnetic catalysis. (The 
talk is based on a series of recent papers with V. Miransky and I. Shovkovy.) This effect 
^ has been established as a universal phenomenon in 2 + 1 and 3 + 1 dimensions: a constant 
^ • magnetic field leads to the generation of a fermion dynamical mass at the weakest attractive 
interaction between fermions 0, |^, ^, 0. The essence of this effect is that electrons behave 
^ i effectively as (l+l)-dimensional ones when their energy is much less than the Landau gap 
\J\eB\ [B is the magnitude of the magnetic field). The lowest Landau level (LLL) plays here 
the role similar to that of the Fermi surface in the BCS theory of superconductivity, leading 
to the dimensional reduction D ^ D — 2 m the dynamics of fermion pairing in a magnetic 
field and to the formation of a chiral condensate at weak coupling. The effect may have 
interesting applications in cosmology p|, |^ and in condensed matter physics 0, as will be 
discussed below. 

The effect of magnetic catalysis was studied in Nambu-Jona-Lasino (NJL) models in 2+1 
0, H and 3+1 dimensions [0, ^, it was extended to the case of external non-abelian chromo- 
magnetic fields |jTO|, finite temperatures [ill] ^^"^ chemical potential [l^, |13], curved space- 



time [0, confirming the universality of the phenomenon. In particular, this phenomenon 
was considered in (3 + l)-dimensional QED 0, ^, |15|, |16|, |1^ |18|. 



We emphasize that we will consider the conventional, weak coupling, phase of QED since 
the dynamics of the LLL is long-range (infrared), and the QED coupling constant is weak 



1 



in the infrared region, therefore, the treatment of the nonperturbative dynamics is rehable 
there. Note that chiral symmetry breaking is not manifested in the weak couphng phase of 
QED in the absense of a magnetic field, even if it is treated nonperturbatively |19]. We will 
show that a constant magnetic field B changes the situation drastically, namely, it leads to 
dynamical chiral symmetry breaking in QED for any arbitrary weak interaction. 
The Lagrangian density of QED in a magnetic field is 



1 

+ 2 



where the covariant derivative is 



D, = d,- leiAl^' + A,), = (o, o) 



(2) 



i.e. we use the so-called symmetric gauge for A^^ . Besides the Dirac index, the fermion 
field carries an additional flavor index a = 1,2, . . . , N. When the bare mass mo = the 
Lagrangian density is invariant under the chiral SUl{N) x SUptiN) x Uv{l) symmetry 
(we will not discuss the dynamics related to the anomalous symmetry Ua{1))- 

We consider first the problem of free relativistic fermions in a magnetic field in 3 + 1 
dimensions and compare it with the same problem in 2 + 1 dimensions. We will see that the 
roots of the fact that a magnetic field is a catalyst of chiral symmetry breaking are actually 
in this dynamics. 

The energy spectrum of fermions in a constant magnetic field is: 



Enip3) = ±^/ml + 2\eB\n + pl n = 0,1,2... (3) 

(the Landau levels). Each Landau level is infinitely degenerate. As the fermion mass rriQ 
goes to zero, there is no energy gap between the vacuum and the lowest Landau level with 
n = 0. The density of states of fermions on the energy surface with Eq = is given by 

\eB\N , , \eB\N 



3 + 1 dimensions, and 



2 + 1 dimensions. 



(4) 



47r2 ' ' ' ~ ' 2tx 

The dynamics of the LLL plays the crucial role in catalyzing spontaneous breaking of chiral 
symmetry. In particular, the density plays the same role here as the density of states 
on the Fermi surface vp in the theory of superconductivity. The next important point is 
that the dynamics of the LLL is essentially (1 + l)-dimensional. Indeed, let us consider the 
fermion propagator in a magnetic field which was calculated by Schwinger 
has the following form (in the chosen gauge): 



20[| long ago and 



S{x, y) = exp ( |(x - yfA^^^x + y) ) S{x - y). 



(5) 



where the Fourier transform of S is 



S{p) 



(pV -pV + ^o)(l + 7^7^tan(eSs)) - pA_-iA_{l + i{m\eBs)) 



(isexp 



^3, ,3 



tan(ei?s) 
eBs 



— m, 







(6) 
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Here p± = (^1,^2), 7± = (71,72) (to get an expression in 2 + 1 dimensions we should put 
Ps = in (D). The propagator S{p) can be decomposed over the Landau level poles as 
follows : 



with 



S{p) = i exp 



Dn{eB,p) = (p|j + mo) 



\eB\ 



E(-i; 



n=0 



Po ~ P3 ~ ^0 ~ '2\eB\n 



(7) 



1 - z7'7'sign(e5))L„ 2 




- (l+^7Vgig^(eS))L„^i(^2^^ +4p^^^Ll_,(^2^^, 

where Ln{x) are the generalized Laguerre polynomials (L„ = L'^-^^{x) = 0). Eq.(|^ implies 
that at rriQ <^ ^\eB\, the LLL with n = dominates and we can write 



S{p) ^ 2i exp( 



pI 



\eB\ 



Pl_+ mo 

Pi 



(8) 



where p|| = p°7°— p'^7^ and = (p'')^ — (p^)^. The matrix O^"-* = (1— z7^7^sign(ei?))/2 is the 
projection operator on the fermion states with the spin polarized along the magnetic field. 
This point and Eq. (H) clearly demonstrate the (H-l)-dimensional character of the dynamics 
of fermions in the LLL. This property is preserved also in the case when the fermion mass is 
generated dynamically. Since at mQ,p|,p5_ ^ \^B\ the LLL pole dominates in the fermion 
propagator, one concludes that the dimensional reduction {D ^ D — 2) takes place for the 
infrared dynamics in a strong {\eB\ >> ttIq) magnetic field. Such a dimensional reduction 
refiects the fact that the motion of charged particles is restricted in directions perpendicular 
to the magnetic field. 

Let us first calculate the chiral condensate in 2 + 1 dimensions for free four-component 
fermions: 



(0|^^|0) = -\imtiS{x,y) = 



AmnN 



"W / d^P [ (is exp 



{2tt)^ J - Jl/A 
\eB\N r°° ds 



tr / d'^pSEip) 
2 2-2 tanh(ei?s) 



-s mg + P3 + p^ 



-mo- 



27r3/2 Ji/A 



r -^e-'^ coth {\eBs\) 



eBs 
eB\N 



271 



mo 0, 



(9) 



where A is an ultraviolet cutoff in Euclidean space and, for concretness, we consider mo > 0. 
Thus, as mo —>■ 0, the condensate {0\i/jijj\0) remains non-zero due to the LLL. Note that the 
expression (^) is nothing else as the Banks-Casher formula relating the fermion condensate 
to the level density of the Dirac operator at zero eigenvalue ||22[. The appearance of the 
condensate in the chiral (flavor) limit, mo — >■ 0, signals the spontaneous breakdown of the 
chiral (flavor) symmetry even for free fermions in a magnetic fleld at D = 2 + 1 [Q]. 
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Repeating the same calculation of the chiral condensate in 3 + 1 dimensions we would 

get 

(0|^^|0) ~ -^o^J^ f In ^ + 0(1) j , mo ^ 0, (10) 

i.e. the condensate is zero and there is no chiral symmetry breaking. Note, however, the 
appearance of logarithmic singularity in ( pil| ) due to the LLL dynamics. As we will see 
below, switching on even a weak attraction between fermions leads to the formation of chiral 
condensate in (3 + l)-dimensional case. 

The above consideration suggests that there is a universal mechanism for enhancing the 
generation of fermion masses by a strong magnetic field: the fermion pairing takes place 
essentially for fermions at the LLL and this pairing dynamics is (1 + l)-dimensional in the 
infrared region. This is the main reason why in a magnetic field spontaneous chiral symmetry 
breaking takes place even at the weakest attractive interaction between fermions in 3 + 1 
dimensions 0, |^. 

Now we shall consider QED in 3 + 1 dimensions whose Lagrangian is given by Eq. (|1|). 
To study chiral symmetry breaking one has to solve the Schwinger-Dyson (SD) equation for 
the dynamical fermion mass. The SD equation for the fermion propagator G{x, y) in an 
external field has the form 

G-\x,y) = S-\x,y) + J:{x,y), (11) 
E{x, y) = AiraY [ G{x, z)T''{z, y, z)V^^{z, x)d'^zd'^z' . (12) 



Here S{x,y) is the bare fermion propagator (^ in the external field S(x, y) is the 

fermion mass operator, and V^yi^x, y), T'^{x, y, z) are the full photon propagator and the full 
amputated vertex. The full photon propagator satisfies the equations 

V-l{x,y) = D-l{x-y) + Ii^,{x,y), (13) 
Iifj,u{,x,y) = -Anati-ff, J d'^ud'^zG{x,u)T^{u, z,y)G{z,x), (14) 

where Dfj,i^{x — y) is the free photon propagator and n^jy(x, y) is the polarization operator. 

It is not difficult to show directly from the SD equations (0), (p^, (0) and ([M]) that 
substitutions 

G{x,y) = exp (^exM^^*(y)) G{x - y), (15) 

r(a;, z) = exp {iex^A'^\y)) f{x - z,y - z), (16) 

T^tJiu{.x,y) =V^^{x -y), (17) 
n^,(x,y) =n^,(x-y) (18) 

lead to equations for translation invariant parts of Green's functions. In other words, in a 
constant magnetic field, the Schwinger phase is universal for Green's functions containing 
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one fermion field, one antifermion field, and any number of photon fields, and the full photon 

propagator is translation invariant. 

We solve the SD equation for the fermion propagator in the so-called ladder approxima- 
tion when the full vertex and full photon propagator are replaced by their bare ones. We 
have 

G{x) = S{x)-Anald''x,d%e'^^^^'^+'^'^^y'^S{x-xi)rG{x,-^^^ (19) 

where the shorthand xA'^^^ly) stands for x^A^^'^ly). 

First, let us show that the solution to the above equation, G{x), allows the factorization 
of the dependence on the parallel and perpendicular coordinates. 

Notice that this form for G{x) is suggested by a similar expression for the bare propagator, 

^(") = 2^^"p(-i)^W^'"^ (21) 



with 

(22) 

^ J (27r)2 fc}-m2' ^ ^ 

taken in the LLL approximation (here / = |ei^|~^/2 is the magnetic length). Performing the 
integrations we arrive at 



g (^xy) = s (^xy) + Awa J -^^d'^xld'^yl exp ^- 



\2 



-tq\\[xi-yi) I s[x" -x\ 



2 

X 7f^(4-Z/})7||5(y!)^^.(g||,g±). (23) 

Regarding this equation, it is necessary to emphasize that the "perpendicular" components 
of the 7-matrices are absent in it. Indeed, because of the identity 0'^~^^^^0^~^ — 0, all those 
components are killed by the projection operators coming from the fermion propagators. By 
switching to the momentum space, we obtain 

(Pll) = (Pll) - / (-%^) - ^")^PM^ (^11' • (24) 

The general solution to this equation is given by the ansatz. 
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where A and B are functions of p|. Making use of this as well as of the explicit form of the 
photon propagator in the Feynman gauge, we get that the function A = 1 while for the mass 
function we get the following integral equation 



mo + 



a 

2^ 



(fqB{q^) f dxexp{-xP/2) 



(26) 



(henceforth we will omit the symbol || from p and q). Thus the SD equation has been reduced 
to a two-dimensional integral equation. Of course, this fact reflects the two-dimensional 
character of the dynamics of electrons from LLL. 

Analytical and numerical analysis of this equation were performed in for the case 
m-o = and in ||18[ for nonzero bare mass. The numerical analysis showed that the so 
called linearized approximation, with B{q^) replaced by the total mass mtot = -B(O) in the 
denominator of Eq. (E^), is an excellent approximation. Then we get 



B{p' 



mo + 



a 



(PqB{q^) 7dxexp{-xiy2) 



27T^ J q^ + m 



tot 



X + (g — pY" 



(27) 



As was shown in 0] (see Appendix C), in the case of weak coupling a and for mo = 0, 
the function B{p) remains almost constant in the range of momenta < p'^<l/P and decays 
like outside that region. To get an estimate for rridyn = -B(O) at a << 1, we set the 
external momentum to be zero and notice that the main contribution in the integral on the 
right hand side of Eq.(p7D is formed in the infrared region with q^<l/P. The latter validates 
in its turn the substitution B{q) B{0) in the integrand of (|26|) , and we finally come to 
the following gap equation: 



5(0) 



a 

2^ 



B{0) 



<Pq 



(ixexp(— Z^x/2) 
q^ + X 



which gives the expression for the dynamical fermion mass (energy gap in the fermion spec- 
trum) : 



rridyn - CVeBexp 



(29) 



where the constant C is of order one and a = e'^/An is the renormalized coupling constant 
related to the scale \/eB. The exponential factor in mdyn displays the nonperturbative nature 
of this result. It can be shown also that the expression (^) for the dynamical mass is gauge 
invariant 0]. 

A more accurate analysis [^, which takes into account the momentum dependence of the 
mass function, leads to the result 



mrf. 



yn 



C^\eB\ 



exp 



TT 



2 V 2a 



(30) 
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The ratio of the powers of this exponent and that in Eq.(|29|) is 7r/2-\/2 ~ 1.1, thus the 
approximation used above is rather rehable. 

We note that mdyn has rather unusual l/\/a behavior of the exponents in (^) and (PI1|). 
Similar dependence was found recently in QCD for a quark gap arising at high densities 
(color superconductivity) |^ . The reason for such a behavior in both cases is the same: the 
presence of long-range interactions. 

To study chiral symmetry breaking in an external field at nonzero temperature we use 
the imaginary time formalism. Now the analogue of the equation (^) (with mo = and the 
replacement mdyn m?{T) in the denominator) reads 



B{u;n',p) = -T Y: / 



dkB{uJn, k) 



dx exp{—Px/2) 



cj2 + p + m2(T)y {iOn - iOn')^ + (k - p)^ + x'- 



(31) 



where Un = vrT(2n + 1) are Matsubara frequencies. 

If we now take n' = 0, p = in the left hand side of Eq.(0) and put B{u!n, k) ^ B{uo, 0) 
const in the integrand, we come to the equation 



dk 



ool + k^ + m^{T) 



dxexp{—Px/2) 

{Un - UJo)^ + k"^ + X 



(32) 



The equation for the critical temperature is obtained putting m(Tc) = and this determines 
the critical temperature |n| 



eB I exp 



mdyn{T = 0), 



(33) 



where rridyn is given by (^). The relationship Tc ~ mdyn between the critical temperature 
and the zero temperature fermion mass was obtained also in NJL model in (2+1)- and (3+1)- 
dimensions (0 and second paper in Ref.[|l^]). The constant C, in the relation Tc = Cmdyn, 
is of order one and can be calculated numerically. We note that the photon thermal mass, 
which is of the order of ^/aT [^, cannot change our result for the critical temperature. 

Taking into account the non-zero bare electron mass we come to the equation for the 
total mass m: 

■log— ^5- ) = mo. (34) 



m cos 




m^ 



It can be shown [T^ that the itterative solution of last equation reproduces all leading double 
logarithmic terms in perturbation theory: 



m = mo 



a , 2 \eB\ 



1 + J- log 



ml 



+ 



a , 2 l^-^l 



24 \ 27r 



log^ 



ml 



+ 



61 



a , 2 1 I 



720 \ 27r 



log^ 



ml 



+ 



(35) 



From Eq. (|3^) we can estimate the dynamical mass due to a magnetic field. For fields of the 
order of ~ lO^^'G which are realized on surfaces of neutron stars we get (m — mo)/m ~ 10%. 
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In the real QED the expansion parameter r] = -^log^deBl/ml) ~ 1 in (|35| ) explores the 
transition between the perturbative regime t] <^ 1 and the nonperturbative massless QED 
regime 1] ^ 1. The value of the parameter r/ ~ 1 is reached at fields of the order ~ lO^^G. 
We recall that strong magnetic fields {B ~ lO^^G) might have been generated during the 
electroweak phase transition It has been speculated in Refs. Q that the character 

of electroweak phase transition could be affected by a generation of a dynamical electron 
mass under such strong fields. The nonperturbative regime becomes prevailing over the 
perturbative one for values of 1] of the order of 2.35 |0 what corresponds to magnetic fields 
~ lO^^G. Ambj0rn and Olesen [^] have claimed that even larger fields, ~ lO^^G, would be 
necessary at early stages of the Universe to explain the observed large-scale galactic magnetic 
fields. 

Since the induced fermion dynamical mass contains an exponential factor (see (p9D, ( |30D ) 
it is quite small at all reasonable values of the coupling a, therefore, there are tiny chances to 
find implications of the magnetic catalysis phenomenon in real experiments. However, it was 
shown recently that the Yukawa coupling and scalar-scalar interaction can considerably 
enhance the fermion dynamical mass (according to the dynamical mass is estimated 



to be rriiiyn — 0.6y\eB\). The most immediate physical implication would be then in the 
electroweak theory. 

Another interesting application of the magnetic catalysis phenomenon is found in (2 + 
l)-dimensional condensed matter systems 0, given the suggestions that high-temperature 
superconductors can be described effectively by (2 + 1) relativistic field theories like NJL 
or QED (the relativistic (Dirac) nature of the fermion fields is related to the fact that 
they describe the quasi-particle excitations about the nodes of a (i-wave superconducting 
gap). According to recent experiments [^, at temperatures significantly lower than T^. of 
superconductivity, the thermal conductivity, as a function of a magnetic field perpendicularly 
applied to the cuprate planes, displays a sharp break in its slope at a transition field 5^, 
followed then by a plateaux region in which it ceases to change with increasing field. The 
critical temperature for appearance of the kink-like behavior scales with the magnetic field 



as ~ ye\B\. This phenomenon may indicate the opening of a second gap, at the nodes 
of the (i-wave superconducting gap, that depends on the strength of the applied magnetic 
field [^, |3^. Indeed, as we saw, in (2 -f- l)-dimensional systems the chiral condensate 
appears even in absence of interaction between fermions. The dynamically generated fermion 



mass scales with a magnetic field like rridyn ~we|S|in2 + l NJL model 0, and rridyn 



alog{y\eB\/a) in QED3 ||28[. The critical temperature for vanishing of the djTiamical mass 
is determined by the dynamical mass at zero temperature (see Eq. (PB|) ) and scales with a 
magnetic field in a way quite similar to the scaling law found in experiments. 

In conclusion, we discuss very breifiy the role of higher order radiative corrections in the 
magnetic catalysis problem. As was shown in Ref.p], because of the (l+l)-dimensional form 
of the fermion propagator of the LLL fermions, there are relevant higher order contributions. 
In particular, considering this problem in the improved rainbow approximation (with the 
bare vertex in the Schwinger- Dyson equations for both the fermion propagator and the 
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polarization operator ), it was shown that the fermion mass is given by Eq. (^) but with 
a a/2. Recently we have shown [^] that there exists a special (non-local) gauge in which 
the SD equations written in the improved rainbow approximation are reliable: in other 
words, in that gauge there exists a consistent truncation of the Schwinger-Dyson equations 
for this non-perturbative problem. The expression for rridyn takes the following form, 

rridyn = C\J\eB\F{a) exp 

where is the number of fermion flavors, F{a) ~ {NaY^^ , Ci ~ 1.82 and the constant C 
is of order one. 

Thus, the magnetic catalysis of chiral symmetry breaking in QED yields a (first, to the 
best of our knowledge) example in which there exists a consistent truncation of the Schwinger- 
Dyson equations in the problem of dynamical symmetry breaking in a (3+l)-dimensional 
gauge theory without fundamental scalar fields. 

I am grateful to the members of the Institute for Theoretical Physics of the University of 
Heidelberg, especially Prof. M.G. Schmidt, for their hospitality during my stay there. This 
research has been supported in part by Deutscher Academischer Austauschdienst (DAAD) 
grant and by the National Science Foundation (USA) under grant No. PHY-9722059. 
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